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A brief review of recent microscopic calculations of nuclear pairing gap is given. A 
semi-microscopic model is suggested in which the ab-initio effective pairing interaction 



is supplemented with a small phenomenological addendum. It involves a parameter 
which is universal for all medium and heavy nuclei. Calculations for several isotopic 
and isotonic chains of semi-magic nuclei confirm the relevance of the model. 
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X 1 1- INTRODUCTION 

M\ 

Up to now, there is no consistent microscopic theory of nuclear matter. The well-known 
Brueckner theory [1] was the first very promising step in this direction but next steps are 
very complicated as one deals with the many-body problem without any small parameter. 
Why the idea to develop the ab initio theory of pairing in finite nuclei is not absolutely un- 
reasonable, although a finite nucleus is much more complicated system than infinite nuclear 
matter? The point is that, for the pairing problem, some simplifications occur in finite nu- 
clei. They originate from the surface nature of nuclear pairing [2] . If the pairing problem is 
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formulated in terms of an effective pairing interaction V^ s in a model space So, this quantity 
turns out to be density dependent [3, 4] with strong dominance of the surface attraction. To 
be more definite, let us write down the simple local 2-parameter ansatz for V£ s within the 
Finite Fermi Systems (FFS) theory [5]: 

V^(ri,r 2 ,r 3 ,r 4 ) = C 

Here C = 300 MeV- fm 3 is the inverse density of state at the Fermi surface, the standard 
FFS theory dimension factor for the effective interaction, and p(r) is the density of the kind of 
nucleons under consideration. Typical values of the parameters (e.g. in [3]) correspond to the 
external constant 7 CX approximately ten times greater than the internal one, 7 111 . Therefore 
it seems reasonable to try to find V^ s starting from the first principles, as the conditions 
for the validity of the Brueckner theory at the surface are much better than inside nuclei. 
Within the Brueckner theory, the gap equation coincides with that of the Bardeen-Cooper- 
Shrieffer (BCS) theory, as the ladder diagrams summation typical of the Brueckner theory is 
made already in the gap equation itself. In any case, the Brueckner theory is valid trivially 
outside the nucleus where all many-body corrections vanish, and therefore it should correctly 
reproduce the 7 ex parameter. This is not the case for the 7 111 parameter (or, more generally, 
for the in-side behavior of the V^(rj) function). But, since it is small, one can hope that 
even noticeable corrections to 7" 1 should not significantly change the gap A value. Such a 
logic has a weak point because of the exponential dependence of the gap on the interaction 
strength which is well known in the weak coupling limit of the BCS theory: 

A F « 25 F exp(l/z/ F \/f ff ) , (2) 

where z/ F = m*k F /-K 2 and e F = fc|/(2m*), m* being the effective mass. This is the reason 
why the knowledge of the in-side behavior of V£ s is important for accurate evaluation of the 
gap value, hence corrections to the BCS theory should be accounted for. 

In the last few years, some progress has been made in the microscopic theory of nuclear 
pairing by the Milan group [6, 7] and Duguet et al. [8, 9]. And some contradictions revealed 
already at the "BCS level", although both the calculations were made within rather close 
frameworks. In particular, the same single-particle spectrum was used for solving the gap 
equation, namely, it was calculated within the Skyrme-Hartree-Fock (SHF) method with 
the Sly4 force which produces the coordinate dependent effective mass m*(r) essentially 



7 ex + (7 in - 7 CX ) 



P(0) 



tf(ri -r 2 )<J(ri -r 3 )£(r 2 -r 4 ). (1) 
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different from the bare one m. In Ref. [7], the value A ~ 1.0 MeV was found for the gap in 
the nucleus 120 Sn (a traditional benchmark for the pairing problem) which is noticeably less 
of the experimental one, A cxp ~ 1.3 MeV. At the same time, in [8] the value A ~ 1.6 MeV 
was obtained for the same nucleus which is essentially larger. In the first case, a lack of 
the gap value is explained in [7] by invoking various many-body corrections to the BCS 
approximation, exchange with low-lying surface vibrations ("phonons") being the main of 
them. Indeed, the latter enlarges the gap value (see, e.g., [6] and [10]) making it closer to 
the experimental value. At the same time, it is rather difficult to find a mechanism that 
can reduce the value of A in Ref. [8]. In Refs. [11, 12] we have analyzed the reasons of 
these contradictions. This point was discussed also in [9]. It turned out that, in fact, these 
two calculations differ in the way they take into account the effective mass. It implies that 
the gap A depends not only on the value of the effective mass at the Fermi surface, as it 
follows from Eq. (2), but also on the behavior of the function m*(k) in a wide momentum 
range. But this quantity is not known sufficiently well [12] that makes rather uncertain the 
predictions of such calculations. To avoid it, we suggest a semi-microscopic model for the 
effective pairing interaction in which the main ab-initio term of V^ s is supplemented with 
a small addendum containing one phenomenological parameter. Preliminary results of this 
model were presented in [13]. 



2. OUTLINE OF THE FORMALISM 

The general form of the many-body theory equation for the pairing gap A reads [5]: 

A = UGG S A, (3) 

where U is the W-interaction block irreducible in the two-paricle channel, and G (G s ) is the 
one-particle Green function without (with) pairing effects taken into account. A symbolic 
multiplication, as usual, denotes the integration over energy and intermediate coordinates 
and summation over spin variables as well. When we used above the term "BCS theory", we 
meant to replace the block U of irreducible interaction diagrams with the free iViV-potential 
V in Eq. (3) and to use the simple quasiparticle Green functions for G and G s (e.g, without 
phonon corrections or others). In this case, Eq. (3) is greatly simplified and can be reduced 
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to the form usual for the Bogolyubov method, 



A = - 



Vx, 



(4) 



where 



J 2ttz 



(5) 



is the anomalous density matrix which can be expressed explicitly in terms of the Bogolyubov 
functions u and v, 



Summation in (6) is carried out over the complete set of Bogolyubov functions with eigen 
energies Ei > 0. 

In Refs. [6, 7], the set of Bogolyubov equations, together with the gap equation (3) 
with the realistic Argonne iViV-interaction v i4 , was solved directly in the basis {A} of states 
restricted to the energy domain up to E max =800 MeV. In addition, as mentioned above, the 
SHF basis with the SLy4 force was used with the coordinate dependent effective mass m*(r), 
which is considerably smaller than the bare mass m. The main difficulty of the direct method 
to solve the nuclear pairing problem comes from rather slow convergence of the sums over 
intermediate states A in the gap equation because of the short-range of the free iViV-force. 
Evidently, this is the reason why the authors of [6, 7] limited the calculations only to one 
nucleus 120 Sn. To avoid the slow convergence problem, the authors of [8, 9] used the super- 
soft "low-k" force Vi ow _k [14] which is defined in such a way that it describes correctly the 
iViV-scattering phase shifts at momenta k<A, where A is a parameter which is not bigger 
than the one corresponding to the limiting energy Eu m ~ 300 MeV, for smaller energy values 
the phase shifts being reproduced accurately. As the force Viow-k vanishes rapidly for k>A, 
one can limit the energy up to £ max ~300 MeV in the gap equation (4). This made it possible 
to calculate in [8] neutron and proton pairing gaps for a lot of nuclei. Usually the low-k force 
is found starting from some realistic AiV-potential V with the help of the Renormaliation 
Group method, and the result doesn't practically depend on the particular choice of V [14]. 
In addition, in [8] Viow-k was found starting from the Argonne potential Vi 8 , which is different 
only a little from the one used in [7], v u . Thus, indeed, the schemes of solving the BCS gap 
equation in [8] and [7] were very similar. 

To overcome the slow convergence problem in the gap equation for finite systems, we used 
a two-step renormalization method. In this approach, we split the complete Hilbert space 
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of the pairing problem S to the model subspace So, including the single-particle states with 
energies less than a fixed value of E , and the subsidiary one, S'. The gap equation is solved 
in the model space: 

A = V^ s GG s A\ So , (7) 

with the effective pairing interaction V£ s instead of the block U in the original gap equation 
(3). It obeys the Bethe-Goldstone type equation in the subsidiary space, 

V^ = U+UGGV^\ S , (8) 

In this equation, the pairing effects could be neglected provided the model space is sufficiently 
large. That is why we replaced the Green function G s for the superfluid system with its 
counterpart G for the normal system. In the BCS approximation, the block U in (8) should 
be replaced by V. To solve equation (8) in non- homogeneous systems, we have found a new 
form of the local approximation, the Local Potential Approximation (LPA). Originally it 
was developed for semi-infinite nuclear matter [15], then for the slab of nuclear matter (see 
review articles [2, 16]) and finally, for finite nuclei [11, 12]. It turned out that, with very 
high accuracy, at each value of the cm. coordinate R, in Eq. (8) the formulae of the infinite 
system embedded into the constant potential well U = U(R) (it explains the term LPA) 
can be used. This simplifies equation for V£ s significantly, in comparison with the initial 
equation for A. As the result, the subspace S' can be chosen as large as necessary. From the 
comparison of the direct solution of Eq. (8) in the slab with the LPA one, it was shown that 
the LPA has high accuracy, even in the surface region, for sufficiently large model space, E 
(~20-=-30 MeV). For finite nuclei (the same 120 Sn), validity of LPA was checked also [11, 12]. 
In this case, the boundary energy should be made larger up to E =A0 MeV. In this article, 
we use the LPA with this value of E for systematic calculations of the gap in spherical 
nuclei. For V, we use just as in [12], the Argonne potential v 18 . 

Let us note that the use of the low-k force VJ ow _k could be also interpreted in terms of 
the two-step renormalization scheme of solving the gap equation (3), with Eq~ 300 MeV 
and with free nucleon Green functions G in (8) (i.e. U(R) = 0). Then, (with U— >V) one 
obtains V£ s — >-Vjo W -k (see [17] where the usual renormalization scheme, similar to ours, is 
used to find Viow-k instead of the Renormalization Group equation). Now, the comparison 
of the direct solution of the gap equation (3) (or (4)) in Ref. [7] with the Argonne NN- 
potential V and of "renormalized" equation (7) with V^ s = VJ ow _k shows that the difference 
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appears because, in the subsidiary subspace S', the effective mass m*^m is used in the 
first case and m* = m, in the second one. Thus, the result for the gap depends not only 
on the value of the effective mass at the Fermi surface, but also on the behavior of the 
function m*(k) in a wide momentum range. This dependence was demonstrated explicitly 
in [11, 12]. The use of the SHF effective mass corresponding to the SLy4 force, or to 
any other version of the Skyrme force, could hardly be approved. Indeed, these effective 
forces were introduced and fitted to describe systematically nuclear masses and radii. As a 
rule, the description of the single-particle spectrum nearby the Fermi surface with Skyrme 
forces is rather poor, and furthermore it is difficult to expect that they will reproduce it 
correctly at those high momenta that are involved in the gap equation (3). This point makes 
it problematic the problem of finding the pairing gap from the first principles completely. 
The situation is even more dramatic because the many-body theory equation (3) contains, 
in addition to the "A;- mass" of the SHF method, the "E-mass" (inverse Z-factor) [18-20], 
which also is not sufficiently well known even in nuclear matter [12]. The corrections to 
the BCS version of Eq. (3) include also the difference of the block U from the polential 
V, mainly due to the so-called induced interaction. The attempt in [7] to find it in terms 
of the same SLy4 force as the nuclear mean field looks questionable. Indeed, this force 
was fitted to the nuclear characteristics which depend mainly on those Skyrme parameters 
determining the scalar Landau-Migdal (LM) amplitudes /, /'. As to the spin amplitudes 
g, g', they remain practically undetermined in the SHF method. At the same time, the 
contribution of the spin channel to the induced interaction is not less than of the scalar 
one [7]. Parameters g, g' are well known from the calculations of nuclear magnetic moments 
within the Finite Fermi Systems (FFS) theory [21], but, just as the Skyrme parameters, 
at the Fermi surface only. But the states distant from the Fermi surface are important to 
calculate the induced interaction. At last, let us imagine to get from some phenomenology 
the functions m*(k),Z(k) and all the LM amplitudes far from the Fermi surface. Even in 
this case, the use of so many phenomenological ingredients devalues significantly the ab initio 
starting point, i.e. the free iViV-potential V in the pairing gap calculation. 

Instead, we suggest to introduce in the effective pairing interaction a small phenomeno- 
logical addendum which embodies, of course approximately, all the corrections to the BCS 
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scheme discussed above. The simplest ansatz for it is similar to Eq. (1) and reads: 

Vefr = K°fF + lC ^5( ri - r 2 )5( ri - r 3 )5(r 2 - r 4 ). (9) 

Here p(r) is the density of nucleons of the kind under consideration, and 7 is a dimensionless 
phenomenological parameter. To avoid any influence of the shell fluctuations in the value of 
p(0), p(0) was averaged over the interval of r<2 fm. The first, ab initio, term in the r.h.s. of 
Eq. (9) is the solution of Eq. (8) with U=V in the framework of the LPA method described 
above, with m*=m in the subspace S'. Then, the gap equation (7) in the model space is 
solved with the self-consistent basis found within the Generalized Energy Density Functional 
(GEDF) method [3] with the functional DF3 where the identity m*=m is assumed. The latter 
is of principal importance for our approach. First, it makes the results less model-dependent, 
all effects of m* 7^ m in both model and subsidiary subspaces being attributed to the in- 
medium corrections beyond the pure BCS approximation. Second, single-particle spectra of 
the GEDF method [3] rule, in better agreement with the experimental ones than 

those of the popular versions of the SHF method [24]. The quality of the single-particle 
spectrum nearby the Fermi surface is very important for obtaining the correct value of the 
gap found from Eq. (4). 

3. ON THE PROCEDURE TO FIND THE "EXPERIMENTAL" GAP 

The gap A is not an observable quantity which can be extracted from experimental data 
directly. Usually, this quantity, A cxp , is found in terms of mass values M of neighboring 
nuclei via 3-term formulae, 

2A+ xp {A) = 5 2 M + = 2M{A + 1) - M(A + 2) - M(A), (10) 

or 

2A cxp (A) = 5 2 M~ = 2M(A - 1) - M(A - 2) - M(A). (11) 
The 5-term expression is usually considered more accurate, being a half-sum of them, 



A cxp (A) = 5 2 M/2 = (5 2 M + + S 2 M~)/2. (12) 

These simple recipes were used, in particular, in [6-9]. However, they originate from the 
simplest model A = const, and the accuracy of such prescription is not obvious a priori. 
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Figure 1. Comparison of the theoretical mass differences 8iM j2 with average gap values Ap for Pb 
isotopes 

To clarify this point we made a calculation which could be considered as a "theoretical 
experiment". We used the GEDF method [3] with the functional DF3 which reproduces the 
mass differences of Eqs. (10), (11) type sufficiently well. First, we calculated the right side 
of Eq. (12) directly, and second, the theoretical gap value. For the latter, we use the "Fermi 
average" combination, 

A f =^(2j+1)A aa /^(2j+1), (13) 

A A 

where the summation is carried out over the states A in the interval of \s\— /i|<3 MeV. A 
similar recipe was used, e.g., in [7]. The comparison of these two quantities is given in fig. 
1 for the lead isotopes and in fig. 2 for the tin isotopes. We see that for the main part 
of nuclei under consideration the difference between values in two neighboring columns is 
within 0.1 MeV. However, there is several cases where it is of the order (or even exceeds) 
0.2 MeV. Leaving aside detailed analysis of these "bad" cases we are forced to put a limit of 
~ 0.1 — 0.2 MeV in the accuracy of the experimental gap determined from Eq. (12). 
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Figure 2. The same as in fig. 1, but for Sn isotopes 

4. CALCULATION RESULTS 



As it was discussed above, for the model space we used the GEDF method by Fayans 
et. al. with the DF3 functional [3]. The model space was extended up to the energy 
Eq = 40 MeV, the subsidiary one up to E max = 1000 MeV. The spherical box of the radius 
R = 16 fm was used, with the grid step h = 0.05 fm. The numerical stability of the results 
was checked by increasing the parameters up to E = 60 MeV, E max = 1200 MeV and 
R = 24 fm, and we found for the gap value a numerical accuracy of 0.01 MeV. 

We calculated the neutron gap for 25 semi-magic isotopes of the lead, tin and calcium 
chains and the proton gap in 9 nuclei, also semi-magic, isotones of the N = 82 chain. The 
formulae above correspond to so-called "developed pairing" approximation [5], i.e. imposing 
the equality of the A + and A~ operators. Therefore we limit ourselves to nuclei having, as 
a minimum, four particles (holes) above (below) the magic core. Therefore, the only isotope 
44 Ca was considered in the calcium chain. 

Let us begin with the neutron pairing. The results are presented in table 1 and figs. 3,4. 
The Fermi average gap values, Eq. (13), found for different values of the parameter 7 in 
Eq. (9). We see that the gap values with the "ab initio" interaction (7=0) are greater by 
30 - 40% than the experimental ones. This difference exceeds significantly the accuracy of 
~ 0.1 — 0.2 MeV for the gap value which we could expect in accordance with discussion of the 
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Figure 3. Neutron gap in Pb isotopes 



> 
CD 



LL 



1.8 
1.6 
1.4 
1.2 
1.0 
0.8 



-•--y = o 
-■- y = 0.07 
a exp 



Sn isotopes 



_l I I I I I I I I I I I I I I I I I I I I I I l_ 



105 110 115 120 125 130 

A 



Figure 4. Neutron gap in Sn isotopes 

previous section. As it can be seen, with few exceptions, it is obtained for 7 = 0.06 — 0.08. 
For the "optimal" value of 7 = 0.07 (the results are exactly half-sums of the values in the 
third and forth columns), the theoretical error exceeds this limit only in 106 Sn and 116 Sn. 
Evidently, it is caused by the fact that the DF3 functional provides an incorrect reproduction 
of the "intruder"-state lh u / 2 , which plays an essential role in the gap equation (7) for these 
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Figure 5. Fermi-average effective pairing interaction 



nuclei. Fig. 3 and fig. 4 are drawn just to illustrate the optimal value of 7. To show that 
the phenomenological addendum to the effective pairing interaction in (9) is indeed rather 
small for 7 = 0.07, we displayed in fig. 5 the localized "Fermi average" of the effective 
interaction. In the mixed coordinate-momentum representation, it is defined as follows: 
V eff (k 1 , k 2 , ri, r 2 ) ->■ V c F ff (i? = r 1 )5(r 1 - r 2 )5(r 1 - r 3 )5(r 2 - r 4 ), where 



V! S (R) = J d 3 tV eS (h = k 2 = k F (R), R - t/2, R + t/2), 



(14) 



with k F (R) = ^/2m(/i — U(R)), provided \x — U(R) > 0, and k F (R) = otherwise. Here /x 
and U(R) are the chemical potential and the potential well of the kind of nucleons under 
consideration. A similar quantity was considered before in the slab system to visualize the 
effective interaction properties [2, 25]. At a glance, the difference between the interaction 
strengths for 7=0 and 7=0.07 is negligible, but it produces noticeable effects in the gap due 
to the exponential behavior in Eq. (2). 

Let us now turn to protons. In this case, the Coulomb potential Vc must be added to the 
expression (9), 

V cff = V: s + V C . (15) 

Again this addendum is small and again it turned out to be important for the gap equation 
due to the the enhancement discussed above. In particular, this was demonstrated in previous 
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Table 2. Proton gap A£ (MeV) for the isotone gap N — 82. 
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calculations [8, 9]. The estimates show that the Coulomb potential could be taken in the 
bare form. Indeed, in the momentum space one has Vc = e 2 /q 2 with a strong maximum at 
small q values provided they persist in the matrix elements < AiA 2 |Vc|A 3 A4 >, with obvious 
notation. In the gap equation, the diagonal elements with Ai = A2 = A3 = A4 are of primary 
importance for which the region around q ~ in the integral dominates. But at small q the 
Coulomb potential Vc{q) is not modified due to the Ward identity. In non-diagonal matrix 
elements the contribution of q ~ kp dominates and Vc{q) could be modified, but in this case 
the contribution of Vc is very small and can be neglected. 

The results for the isotone chain N=82 are given in Table 2 and displayed in fig. 6. To 
demonstrate the effect of the Coulomb interaction, we show the results with the interaction 
Vg ff = Vg ff which difference from the corresponding value in column 3 gives exactly the 
Coulomb effect in the gap. Indeed, it is rather big (about 0.5 MeV), in qualitative agreement 
with [8]. Again at 7=0.07 the agreement is almost perfect for the most part of nuclei, 
and only for the two heaviest isotones the disagreement exceeds 0.2 MeV. In this case, 
the possible reason lies in the proximity to the phase transition to the deformed state (at 
A ~ 150). Average difference between the theoretical and experimental gap values for 34 
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Figure 6. Proton gap in N = 82 isotones 



nuclei considered is equal to \J (t)A) 2 ~0.13 MeV. As it follows from the analysis in Sect. 
3, this value is within the accuracy of the experimental values of the gap defined with the 
relation (12). 



5. CONCLUSIONS 



We suggest a simple semi-microscopic model (9) for the effective pairing interaction 
containing one phenomenological parameter which takes into account approximately var- 
ious corrections to the pure BCS theory. This model reproduces rather well experimen- 
tal values of the neutron and proton gaps in semi-magic nuclei. The overall agreement 
{\J (5A) 2 ~0.13 MeV) is better than that obtained in [8], where the authors did not intro- 
duce free parameters explicitly but they made it implicitly by using a specific k-dependence 
of the effective mass. 

The ansatz of Eq. (9) possesses an obvious drawback. The phenomenological GEDF 
pairing interaction of [3] contains the surface term ((x(dp/dr) 2 ) that plays an essential role 
for the description of the odd-even effect (staggering) in nuclear radii. It originates mainly 
from the exchange by surface phonons which was explicitly taken into account in [6, 7]. The 
addition of such a term in Eq. (9) is associated with introducing a new parameter, and at 
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the first stage we preferred to avoid it. A more consistent scheme should, evidently, include 
the explicit consideration of the low-lying phonons, as e.g. in [6], but taking into account 
the so-called tadpole diagrams [23]. In this case, the phenomenological constant 7, of course, 
will change. 
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